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Q ■ Abstract 

Q I The introduction of a chirally twisted mass term has been proposed as an attractive 

I approach to 0(a) improvement of Quantum Chromodynamics with Wilson fermions on a 

I lattice. For numerical simulation projects it is important to know the phase structure of 

rin ! the theory in the region of small quark masses. We study this question in the framework 

Qh! of chiral perturbation theory. Generalizing the analysis of Sharpe and Singleton we find 

I extensions of normal and Aoki phase scenarios and a possible new phase with spontaneous 

^ I breakdown of chiral symmetry to a discrete Z2. 
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1 Introduction 

In numerical simulations of Quantum Chromodynamics on a lattice one of the most important 
issues is to reduce lattice artifacts in order to improve the approach to the continuum limit. 
Recently the lattice formulation of QCD with Wilson fermions and a chirally twisted quark 
mass matrix E] has attracted much interest in this respect. Frezzotti and Rossi |31 S] have 
advocated this approach as a way to remove lattice artifacts to first order in the lattice spacing a 
by a suitable choice of the chiral twist. Exploratory numerical studies [HI have shown favourable 
algorithmic features of simulations of twisted mass lattice QCD. 
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A prerequisite to any numerical simulation project is the knowledge of the phase structure 
of the model under consideration. Where are lines or points of phase transitions of first or 
higher order and how do physical quantities like particle masses behave near them? It is the 
purpose of this note to discuss the phase structure of twisted mass lattice QCD on the basis 
of chiral perturbation theory. Similar considerations have been made independently by Sharpe 
and Wu P[Z|. 

We consider lattice QCD with Nf = 2 quark flavours and degenerate quark masses m. Let 
rric be the "critical" value of m, where the quark condensate has a discontinuity, and 

m = m — rric (1) 

be the subtracted quark mass. The chirally twisted mass term can be introduced in the socalled 
"twisted basis" in the form g(mc + M{u!))q with 

M{uj) = rriqe'^^'^' =m + i/j-f^n , (2) 

where 

rh = rriq cos(ci;), fj, = rrig sin(a;) , (3) 

and Ta denote Pauli matrices. 

The low-energy regime of QCD is described by chiral perturbation theory [HI El [El in terms 
of the dynamics of the pseudo-Goldstone bosons of broken chiral symmetry, which are the pions 
in our case. The pion fields 7r{,(x) enter the SU(2)-valued matrices 

U{x) = exp '^bix)Tb^ , (4) 

which are the fundamental variables of chiral perturbation theory. The above representation in 
terms of pion fields refers to the vacuum U = 1 and is modified accordingly if the vacuum is 
at a different point. In leading order the chiral effective Lagrangian including lattice effects up 
to 0{a) is dllini 

£2 = ^ Tr {d,U^d,U) - ^ Tr (^f/^ + Ux^) - ^ Tr (pt/t + f/pt) , (5) 

where 

X = 25o(ml - iprg), (6) 
appears in the chiral symmetry breaking mass term, and 

p = 2Woa 1 (7) 

parameterizes the lattice artifacts. The effective Lagrangian in next to leading order including 
lattice artifacts up to order has been calculated by Bar, Rupak and Shoresh The 
introduction of a twisted mass term in chiral perturbation and its application to the calculation 
of pion masses and decay constants have been made in [T2l fT4l . 

2 



The discussion of the phase diagram of lattice QCD in the framework of chiral perturbation 
theory has been pioneered by Sharpe and Singleton [1^]. It is based on the potential contained 
in the chiral Lagrangian. Let 

U = Uol + iUaTa, a = 1, 2, 3, (8) 

so that 

M = (mO,Mi,U2,M3) (9) 

is a unit 4- vector: u ■ u = 1. The physical significance of the variable Uq is given through its 
relation to the chiral condensate 

{qq) = -2F^Bo{uo) (10) 

in leading order. The potential in next to leading order for the case of vanishing twist is of the 
form 

V = -CiUo + C2ul, (11) 

where 

Cl = 2F2(5om, + Woa) . (12) 

Sharpe and Singleton discuss the region of quark masses, where they are small of order a^, in 
which case both coefficients ci and C2 are to be considered of order a^. 

In order to study the phase diagram in the niq — fi plane the effect of the twisted mass term 
has to be incorporated. Its addition to the QCD Lagrangian amounts to extend the potential 
to 

V = -Cl Uo + C2 Mq + C3 U3 , (13) 

where 

c, = 2F^Bofi. (14) 

As for the quark mass, the twist mass /i is also counted as being of order a^, so that higher 
powers of /i are neglected for the moment. 

Depending on the sign of C2 there are now two possible scenarios for the phase structure, 
which will be discussed in the following. 

2 Aoki scenario 

We start by considering the case C2 > 0, where an Aoki phase [H] appears. Let us first briefiy 
recall the analysis of Sharpe and Singleton [111 for /i = 0. If 

|e| > 1, where e = -r— r, (15) 

2|C2| 

the potential has its minimum outside the allowed region \uq\ < 1 and consequently the vacuum 
is at Mo = 1, ^7 = 1 or at Mo = — 1, = — 1, depending on the sign of ci. The pion mass is 
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given by 



ml = ^{\c,\-2c2). (16) 

Chiral symmetry SU(2)^ ® SU(2)j:j is broken explicitly to flavour SU(2)y. 
On the other hand, for small quark masses, namely in the range 

|e| < 1, (17) 

there is a small region of Aoki phase. The potential has its minimum at = e and the remaining 
components of u develop a vacuum expectation value associated with spontaneous flavour and 
parity breaking, which can be chosen to be directed along the third axis: = ±Vl — e^- The 
nonvanishing Us corresponds to 

ixi^Tsx) ^ 0. (18) 

The remaining flavour symmetry is U(l) and tti and 712 are the massless Goldstone bosons. The 
neutral pion is massive with a mass 

ml,= ^-0{l-e'). (19) 

Varying Ci inside the Aoki phase from e=ltoe = — 1, the chiral condensate changes continu- 
ously from negative to positive values. 

The boundary of the Aoki phase is at |e| = 1, where there is a second order phase transition 
and all pions are massless. The width of the phase is proportional to a^. 

Now let us switch on the twist mass to a nonvanishing value, which we may choose to be 
positive: /x > 0. This leads to a shift of the minimum of the potential in the third direction. As 
a consequence, for all quark masses, M3 acquires a nonvanishing negative value, which is given 
by the solution of a quartic equation and is shown in Fig. [H 

Instead of writing down the exact solution in terms of the formulae of Ludovico Ferrari ^\ 
here, we restrict ourselves to displaying it to first order in C3 oc /i for |e| away from 1: 



+ Oicl) for |e| > 1 



|ci| - 2o 



«3= ' " ,2 ^ (20) 

-v/r^-^^ + 0(c^) for |e|<l. 

The breaking of flavour SU(2) is in this case explicit, and in addition to 713 also the charged 
pions TTi and n2 are massive. 

Uq changes continuously as a function of the quark mass and there is no phase transition. 
The phase transition, that existed for /i = 0, is washed out for /i 7^ 0. The phase diagram is 
shown in Fig. [21 Between the second order endpoints is a line of first order transitions, where 
U3 makes a jump when going from positive to negative values of /x. 
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Figure 1: The values of Uq and M3 in the minimum of the potential as a function of Ci/2|c2| oc m 
for C3 = C2 in the Aoki scenario. 
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Figure 2: Phase diagram for the Aoki scenario. 

3 Normal scenario 

A different scenario results for C2 < 0. The case of vanishing /i has been discussed by Sharpe 
and Singleton jTHj too. For fi = the potential has its minimum at = 1 for Ci > 0, i.e. for 
positive quark mass, and at Uq = —1 for ci < 0. At Ci = 0, corresponding to the quark mass 
being equal to rric, there is a first order phase transition, where Uq jumps from +1 to —1. The 
chiral condensate makes a jump of size —AF^Bq. The pion masses are nonzero and are given 

by 

ml = ^{\c,\+2\c2\). (21) 

For /i 7^ the minimum of the potential is shifted away from Uq = ±1 by the "magnetic field" 
/i. We have |mo| < 1 and M3 7^ 0. The corresponding fiavour and parity symmetry breaking is 
explicit and the pions remain massive. 

The values of uo and U3 are again given by the solution of a quartic equation, and are shown 
in Fig. O 
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Figure 3: The values of Uq and M3 in the minimum of the potential as a function of Ci/2|c2| oc m 
for C3 = \c2\ in the normal scenario. 



To first order in /i we have 

U3 



C3 



|ci| + 2|C2 



+ Oicl) 



2(|ci|+2|c2|)2 



0{cl). 



(22) 
(23) 



At Cl = there is a first order phase transition, where uq changes sign discontinuously, as 
indicated in Fig. HI The value of M3 at the transition is 



^^3 



C3 



2|C2| 



(24) 



For increasing /z the minima of the potential at the transition point Ci = approach each 
other. There is an endpoint of the transition line at 



C2 2 

~ a 



(25) 



and the phase diagram looks like it is shown in Fig. O 

On the phase transition line the neutral pion mass behaves as 

1 



■{44 - 4) ~ a' 



(26) 
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Figure 4: Jump of Uq at the phase transition for the normal scenario. 
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Figure 5: Phase diagram for the normal scenario, 
and the same holds for the jump of the condensate, 



1 - 



4_ 

4ci 



2\C2\ 



m^3 , 



(27) 



as displayed in Fig. El 

The masses of the charged pions are continuous across the first order line, where they have 
their minimum 

2|C2| 2 



mil = 



~ a 



(28) 



4 A new phase 

So far we have considered the case of twist masses /i of order a? with only a linear term in M3 in 
the potential. In general, higher terms in the Sharpe-Singleton potential have to be included. 




Figure 6: Pion masses and chiral condensate as a function of /i along the phase transition line 
in the normal scenario. 

At the next order additional terms quadratic in Ua have to be added: 

V = -ci uq + C2 Mq + C3 U3 + C4 M3 + C5 U0U3. (29) 

What are their possible consequences? In the normal scenario, C2 < 0, one finds that these 
terms do not lead to qualitative changes of the phase diagram. Only numerical modifications 
will occur. 

On the other hand, in the Aoki case the new terms can lead to new phases, if the coefficient 
C4 oc /i^ is positive. In this situation a minimum of the potential in the interior of the Uq-u^ unit 
circle will occur for sufficiently large fi. As a consequence, fiavour symmetry is broken down 
spontaneously to a discrete Z2 by nonvanishing expectation values of ui or U2- In addition 
to the segment of Aoki phase at /i = 0, regions of this new phase in the upper (/i > 0) or 
lower (yU < 0) half-planes can form. Inside this phase the pion masses vanish and the quark 
condensate varies continuously. 

Along the lines of catastrophe theory [TH], terms in the potential of higher orders than 
quadratic do not change the discussed scenarios for small masses qualitatively as long as none 
of the quadratic coefficients above vanishes accidentally. 

5 Discussion and conclusion 

As we have seen, there can exist a line of first order phase transitions in the mq-fi plane, 
terminating in second order endpoints. Depending on the sign of the coefficient C2, the line is 
horizontal or vertical in the plane. Whereas for values of the inverse gauge coupling /3 < 4.6 the 
Aoki scenario appears to be realized [in], recent Monte Carlo calculations at /3 = 5.2 indicate 
the presence of the normal scenario. The first order phase transition line for nonvanishing /x 
occurs at a twist angle of 7r/2, where 0{a) improvement is predicted [Sill]- The associated two- 
phase coexistence and metastability hampers numerical calculations. Nevertheless, interesting 
physical observables, like decay constants and squared masses, are continuous across the phase 
transition line and are expected to be 0{a) improved there. 
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From the point of view of Monte Carlo calculations it would be desirable to have the phase 
transition line as short as possible. There should be a value of (5 somewhere between 4.6 and 
5.2 where jic oc C2 vanishes. This point, however, would not be close enough to the continuum 
limit. In order to make jic small one could go nearer to the continuum limit at larger values 
of /3, or try to improve the action in such a way that C2 is small. One can search for such a 
situation on the untwisted axis (// = 0) by monitoring the minimum of the pion masses and 
trying to make it small. The jump of the condensate is not suitable for the search, as it is 
insensitive to C2. 
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References 

[1] R. Frezzotti, RA. Grassi, S. Sint and R Weisz, Nucl. Phys. B (Proc. Suppl.) 83 (2000) 
941; JHEP 0108 (2001) 058. 



[2] 


R. 


Frezzotti, Nucl. Phys. B (Proc. Suppl.) 119 (2003) 140. 


[3] 


R. 


Frezzotti and G.C. Rossi, hep-lat/0306014. 


[4] 


R. 


Frezzotti and G.C. Rossi, hep-lat/0311008. 


[5] 


F. 


Farchioni et al., hep-lat/0406039, and private communications. 


[6] 


S. 


Sharpe and J. Wu, private communications. 


[7] 


J. 


Wu, talk at the conference LATTICE2004, June 23, 2004, and article in preparation 


[8] 


s. 


Weinberg, Physica A 96 (1979) 327. 


[9] 


J. 


Gasser and H. Leutwyler, Ann. Phys. 158 (1984) 142. 


[10] 


J. 


Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 


[11] 


G 


Rupak and N. Shoresh, Phys. Rev. D 66 (2002) 054503. 


[12] 


G 


Miinster and C. Schmidt, Europhys. Lett. 66 (2004) 652. 


[13] 





Bar, G. Rupak and N. Shoresh, hep-lat/0306021. 


[14] 


G 


Miinster, C. Schmidt and E. Scholz, hep-lat/0402003. 


[15] 


S. 


Sharpe and R. Singleton, Phys. Rev. D 58 (1998) 074501. 



9 



[16] S. Aoki, Phys. Rev. D 30 (1984) 2653; Phys. Rev. Lett. 57 (1986) 3136. 

[17] in: G. Cardano, Artis magnae sive de regulis algebraicis, Liber unus, Joh. Petreium, Niirn- 
berg, 1545. 

[18] A. Held and P. Yodzis, Gen. Rel. and Grav. 13 (1981) 873. 

[19] E.M. Ilgenfritz, W. Kerler, M. Miiller-Preuftker, A. Sternbeck and H. Stiiben, Phys. Rev. 
D69 (2004) 074511. 



10 



